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Abstract 



This technical report is divided into two parts. The first part of the technical report presents a novel framework 
for modeling the uplink and downlink intercell interference (ICI) in a multiuser cellular network. The proposed 
framework assists in quantifying the impact of various fading channel models and multiuser scheduling schemes 
on the uplink and downlink ICI. Firstly, we derive a semi-analytical expression for the distribution of the location 
of the scheduled user in a given cell considering a wide range of scheduling schemes. Based on this, we derive the 
distribution and moment generating function (MGF) of the ICI considering a single interfering cell. Consequently, 
we determine the MGF of the cumulative ICI observed from all interfering cells and derive explicit MGF expres- 
sions for three typical fading models. Finally, we utilize the obtained expressions to evaluate important network 
performance metrics such as the outage probability, ergodic capacity and average fairness numerically. Monte-Carlo 
simulation results are provided to demonstrate the efficacy of the derived analytical expression|3- 

The second part of the technical report deals with the statistical modeling of uplink inter-cell interference (ICI) 
considering greedy scheduling with power adaptation based on channel conditions. The derived model is utilized 
to evaluate important network performance metrics such as ergodic capacity, average fairness and average power 
preservation numerically. In parallel to the literature, we have shown that greedy scheduling with power adaptation 
reduces the ICI, average power consumption of users, and enhances the average fairness among users, compared 
to the case without power adaptation. 



Explosive growth in the demand of high quality wireless data services compel the network designers 
utilize spectrum more aggressively which on one side enhances the spectrum efficiency, whereas 

'The first part of the technical report is currently submitted to IEEE Transactions on Wireless Communications 
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on the other side it enhances the intercell interference (ICI) which is an alarming bottleneck in the 
telecommunication growth paradigm. The allocation of the same frequency bands across neighboring 
cells produces indeterministic ICI which is highly dependent on the statistics of the channel characteristics 
and on the dynamics of the multiuser scheduling decisions. In this context, it is of immense importance 
for the system designers to accurately characterize and investigate the behavior of the ICI which helps 
in gaining more theoretical insights, quantifying various network performance metrics and developing 
efficient resource allocation and interference mitigation schemes. 

Orthogonal frequency division multiple access (OFDM A) has been recently adopted as the multiple 
access scheme for the state-of-the-art LTE and WiMAX cellular technologies. In OFDMA, a wide-band 
frequency- selective fading channel is decomposed into a set of orthogonal narrow-band subcarriers. The 
orthogonality among the subcarriers per cell makes the intra-cell interference almost negligible. However, 
with universal frequency reuse among cells (i.e., all cells use the same set of subcarriers), the ICI at each 
subcarrier may cause severe degradation in the network performance. In OFDMA networks, the subcarriers 
are allocated adaptively among users per cell based on a predefined scheduling scheme. Moreover, each 
subcarrier is allocated to only one user per cell at a given time instant and, thus, the number of interfering 
users on each subcarrier is rather limited. Therefore, the cumulative ICI on a given subcarrier may not be 
modeled accurately as a Gaussian random variable (RV) by simply invoking the central limit theorem. 

Several recent studies considered the modeling of ICI in the downlink where the location of interferers 
is typically deterministic. In [|T|, a semi-analytical approach to estimate the ICI by modifying the Burr 
distribution considering path loss, Rayleigh fading and log-normal shadowing was presented. A semi- 
analytical distribution for the signal-to-interference-noise ratio (SINK) has been derived in [|2l under path 
loss and log normal shadowing for randomly located femtocell networks. In [3], the applicability of the 
Gaussian and binomial distributions for modeling the downlink ICI is investigated. In [4j, the optimal 
threshold is derived for fractional frequency reuse (FFR) systems assuming ICI as Gaussian RV. In 
the authors derived the distribution of the ICI under log-normal shadowing and Rayleigh fading. The 
distribution of ICI is shown to highly deviate from the Gaussian distribution in OFDMA networks. 

In comparison to downlink, the nature of uplink intercell interference (ICI) is different in various 
aspects that include the following: (i) Due to the implicit symmetry and fixed locations of the BSs in the 
typical grid-based downlink network models, the number of significantly contributing interferers typically 
remains the same irrespective of the position of the mobile receiver. Also, it has been shown in [|3l that 
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the strongest interference is generated by two closest interfering BSs irrespective of the mobile receiver 
location. However, the number of significantly contributing interferers in the uplink cannot be quantified at 
a given instant due to the highly varying locations of the interfering mobile transmitters; (ii) Conditioned 
on the location of the desired mobile receiver within a cell, the exact distance of the interfering BSs can 
be calculated in the typical grid-based downlink network models. However, knowing the location of the 
BS receiver in the uplink does not help in determining the exact location of the interfering mobile users; 
(iii) In the uplink, cell edge and cell center mobile users are subject to the same amount of interference 
on a given subcarrier, which is the interference received at the BS. Whereas the same is not true for the 
downlink in which cell edge users experience higher interference coming from the nearby BSs flU, flT]. 

From the system design perspective and performance analysis, both uplink and downlink are equally 
important for the network designers. However, based on the above arguments, the nature of the uplink 
ICI on a given subcarrier is different from the downlink ICI due to the inherit differences in the nature 
of the interferers. Therefore, uplink and downlink require different modeling approaches to statistically 
capture their interference dynamics. In this context, the proposed analytical approach to model uplink ICI 
on a given subcarrier is novel as it captures the impact of various state-of-the-art scheduling schemes and 
is generalized for different composite fading models. The approach provides the MGF of the ICI which 
can be utilized to calibrate various system performance metrics. 

Some worth mentioning research works for the uplink appear in H, [|8l- [[T0l . In H, the authors 
developed an analytical model for subcarrier collisions as a function of the cell load and frequency reuse 
pattem. They derived an expression for the SINR in the uplink and downlink, ignoring the effect of 
shadowing and fading. In [8], the authors developed an analytical expression for the subcarrier collision 
probability considering non-coordinated schedulers. In the authors modeled uplink ICI in an OFDM A 
network as a function of the reuse partitioning radius and traffic load assuming arbitrary scheduling. In 
ifTOl . the authors presented a semi-analytical method to approximate the distribution of the uplink ICI 
through numerical simulations without considering the impact of scheduling schemes. 

In this paper, we propose a novel theoretical framework to derive the statistics of the uplink and downlink 
ICI on a given subcarrier as a function of both the channel statistics (i.e., path loss, shadowing and fading) 
and multiuser scheduling decisions. The framework is generic in the sense that the derivations hold for 
generalized fading channels and various scheduling algorithms. We start by deriving the distribution of 
the location of the scheduled user in a given cell. We then derive the distribution and moment generating 
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function (MGF) of the ICI considering a single interfering cell. Next, we derive the MGF expression for 
the cumulative ICI experienced from all interfering cells over generalized fading channels, and present 
explicit expressions for three practical fading models. Finally, we demonstrate the importance of the 
derived expressions by utilizing them to evaluate important network performance metrics. 
Notation: Exp(A) represents an exponential distribution with parameter A, Gamma(ms, mc) represents 
a Gamma distribution with shape parameter and scale parameter uic. /Cg("^c, "^s, fi) represents the 
Generalized- /C distribution with fading parameter uic, shadowing parameter and average power Cl. 
r(.) represents the Gamma function. p{A) denotes the probability of event A. /(.) and F{.) denotes 
the probability distribution function (PDF) and cumulative distribution function (CDF), respectively, [a, b] 
denotes a discrete set of elements which ranges from a to b. Finally, ¥.[.] denotes the expectation operator. 

II. Part I: System Model and Proposed Framework 
A. Description of the System Model 

We consider a given cell surrounded by L interfering neighboring cells. For analytical convenience, the 
cells are assumed to be circular with radius R. Each cell / is assumed to have U uniformly distributed 
users. The frequency reuse factor is assumed to be unity with each subcarrier reused in all cells. The 
bandwidth of a subcarrier is assumed to be less than the channel coherence bandwidth, thus, each subcarrier 
experiences flat fading. Time is divided into time slots of length smaller than the channel coherence time 
and, thus, the channel variation within a given time slot is negligible. 

Generally, the scheduling strategies can be broadly categorized into two classes; (i) rate maximization 
(i.e., rate adaptation) while transmitting with constant/maximum power; (ii) power minimization (i.e., 
power adaptation) while achieving a fixed data rate. In this work, we focus on rate adaptive schemes 
where users transmit with their maximum power in order to maximize their rate depending on the existing 
channel and interference conditions. Therefore, for the scope of this paper, we assume that all users 
transmit with their maximum power Pmax on a given subcarrier with rate adaptation depending on their 
channel qualities. At this point, it is also important to emphasize that this is not a limitation and the 
approach can be extended for various uplink power control mechanisms. The instantaneous signal to 
noise ratio (SNR) 7 of each user can then be written as follows: 

7 = PmaxC^ = Kr-^C (1) 



where K = ^•^^^'^ , C is the path loss constant, r is the user distance from its serving BS, il) and 77 denotes 
the shadowing and small scale fading coefficient between user and BS on a given subcarrier, respectively, 
(3 is the path loss exponent, a denotes the thermal noise at the receiver and C, is the composite fading. 
Note that all users are assumed to be associated with their closest BS (611 , IfTTTl . therefore r < R. 

Each cell is divided into K concentric circular regions. Since path loss decays exponentially from cell 
center to cell edge, therefore, we consider discretization of cellular region in such a way that the path 
loss decay within each circular region remains constant or uniform. The main motivation for dividing the 
cell into a discrete set of circular regions relies on the fact that the channel statistics of the users located 
in a given circular region become relatively similar especially for large values of K. More explicitly, the 
characterization of the circular regions can be demonstrated as follows: 

log,„r. = ^±i5^^. r<R (2) 

where k is the path loss decay within each circular region [dB]. Due to the exponential nature of the 
path loss, it varies rapidly near the cell center than at the cell edge, therefore, (i) each of the k^^ circular 
region bounded by two adjacent rings, i.e., and rk-i possess non-uniform width = — rk-i', (ii) 
the number of circular regions are high in the cell center than at the cell edge; (iii) the average number 
of users located within k^^ circular region bounded by ring and Vk-i are considered to be located at 
Tfc. Note that, this is an approximation which is required for deriving the analytically tractable model of 
ICI and in any case it is not required for the Monte-Carlo simulations. The average number of users in 
each ring k (for analysis) can then be given as: 

U{rl — r^_i) 7-10 
"fc = ^2 k = l,2,--- ,K (3) 

It is important to note that can be a fraction; therefore, we round off the fractional part of users. 

B. Main Steps of the Proposed Framework 

In order to characterize the statistics of the uplink and downlink ICI for generalized fading channels 
and various scheduling schemes, the proposed framework mandates the following steps: 
i) Derive the distribution fr^^^{r) of the distance of the allocated user in a given cell from its 
serving BS based on the deployed scheduling scheme. Without loss of generality, we assume the 
same scheduling scheme is implemented in all cells; therefore, fr^^iir) remains the same for all cells. 
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ii) Derive the distribution (f ) of the distance between the allocated user in a neighboring interfering 
cell and the BS of the cell of interest fsci- 

iii) Derive the distribution fxiix) of the interference from the neighboring cell /. Finally, derive the MGF 
of the cumulative ICI, i.e., Y = J2i=i-^i- 



Considering the high dependence of the uplink ICI on the location of the scheduled users in the neighbor 
interfering cells which in turn depends on the deployed scheduling schemes, we derive in this section 
the distribution of the distance between the scheduled user and its serving BS in a given cell (i.e., the 
probability mass function (PMF) of Tsci) considering the following five scheduling algorithms: greedy 
scheduling, proportional fair scheduling, round robin scheduling, location based round robin scheduling, 
and greedy round robin scheduling. 

A. Greedy Scheduling Scheme 

Greedy scheduling is an opportunistic scheme that aims at maximizing the network throughput by taking 
full advantage of multiuser diversity. However, it suffers from low fairness among users which makes it 
less attractive for network operators. The procedure for determining the PMF of Tsci considering greedy 
scheduling is divided into two steps: 

Step 1 (Selecting the user with the highest SNR in ring k): Since the path loss decay within each circular 
region is considered to be uniform, we approximate the distance of all users located within k^^ circular 
region by ring for analytical tractability as we already mentioned in Section II. In this step, we select 
a user with maximum SNR in each ring k which posses users. Thus, selecting a user in a ring k is 
equivalent to selecting the user with maximum channel gain among all the users in ring k, i.e.. 



where Q is the composite fading channel gain between user i and its BS on a given subcarrier. The CDF 
and PDF of the maximum channel gain Ck can be written as follows: 

F^ACk) = \[FdCky= {F^{Ck)r (5) 

i=l 



III. Part I: Distribution of the Scheduled User Location 



Cfc = max{Ci,C2,-- - ,0, ■■■ ,C«J 



(4) 




(6) 



To consider path loss, we now perform a transformation of RVs using ([T]), 7^ = Kr^^ ^C^, where, 7^ is 
the selected user SNR in each ring k. The CDF and PDF of 7^ can then be written as follows: 

i=l 

fM = ^E4(7.-fi^-^) n F,X7,4k~^) '= —.fdl^rlK-^) (i^c(7.rfi^-^)) '~ (8) 

j=i i=l,i^j ^k 

Step 2 {Selecting the user with maximum SNR among K rings): In this step, we compute the probability 
of selecting the k^^ ring among all other rings. It is important to note that this is equivalent to selecting 
the ring k which possesses the user with the highest SNR among all rings. Conditioning on 7^, the PDF 
of Tsei can be written explicitly as follows: 

K K 

P(rsei = rfc|7fc)= n V{l^<lk)= n ^7.(7fc) (9) 

i=l,ij^k i=l,ij^k 

By averaging over the distribution of 7fc, the final expression for the PMF of r^ei is 

l-oo f K \ 

P(rsei = rk)= n fyMd^k (10) 

•^0 \i=i,i^k J 

Using dV]), (flOl) can be written for i.i.d. case as follows: 



"'O i=i^i^k ^k 

where r^d G [0,R]. The results in (fTTI) are generalized for any shadowing and fading statistics. Even 
though (fTT)) is not a closed form expression, the integration can be solved accurately using standard 
mathematical software packages such as MAPLE and MAT HEMATIC A. 

B. Proportional Fair Scheduling Scheme 

The proportional fair scheduling scheme allocates the subcarrier to the user with the largest normalized 
SNR (7/7) lfT2l . where 7 and 7 denote the instantaneous SNR and the short term average SNR of a given 
user, respectively. In other words, the selection criterion is based on selecting a user who has maximum 
instantaneous SNR relative to its own average SNR. The distribution of rgei can be derived as follows: 
Step 1 (Selecting the user with maximum normalized SNR in ring k): In this step, the performance of 
proportional fair scheduling scheme is independent of the path loss factor if users are moving relatively 



slowly, i.e., their path loss remains nearly the same on a short term basis. In this case, the problem of 
selecting the maximum normalized SNR in a ring k can be written as: 

Cfc = max <^ ^, ^, ■ ■ ■ , ^, ■ • • , ^ ^ (12) 

where = /o°° 0/Ci(0)'^0 is the average of the composite fading channel and Cfc is the maximum 
normalized composite fading channel gain in ring k. For i.i.d. composite fading gains of the users located 
in ring k, i.e., C = Ci = C2 = ■ ■ ■ Cwfe' the problem of selecting the user with maximum normalized channel 
gain reduces to selecting the user with the maximum channel gain, i.e., Ck = max {Qi, C2, ■ ■ ■ , 0; ' ' ' ; Cnfe}- 
Thus, for any ring k, the CDF and PDF of the selected SNR 7^ = Kr~^^C,k can be written as: 

^7.(7.) = \{FdC^lkrlK~') '= {FaClkrtK-')r (13) 

i=l 

Uk Uk 

fM = — E4(07.rfK-^) n FUQlkr'.K-') '= ^{FUClkr',K-')r-'fdClk4K-') 

^k j=l i=l,i^'j 

(14) 

The short term average SNR of the selected user in ring k, i.e., 7^ can then be computed as 'jk = 
Jo° lkf'yki'lk)d'yk- Finally, the normalized selected SNR in each ring k can be defined as = ^ and 
performing a transformation of RVs, the CDF and PDF of can be given as F^^{^k) = -^7^(7^^^); and 
hki^k) = ^J^kilkik), respectively. 

Step 2 {Selecting the ring k with maximum normalized SNR from the K rings): Once we characterize the 
PDF and CDF of S,k, the probability of selecting any ring k can be written using ^ as follows: 

noo f K \ 

P(rsei = rk)= n ^5.(^0 hk{ik)dik (15) 

"^0 \i=i,i^k J 

C. Round Robin Scheduling Scheme 

Round robin scheduling is a non-opportunistic scheme where a user is selected randomly within a time 
slot. As each user has equal probability of allocation, the scheme can be referred to as a strictly fair 
scheduling scheme. The round robin scheme provides maximum fairness among users and may serve as 
a lower bound in terms of network throughput which is useful in calibrating the performance of other 
scheduling schemes, however, the resulting network throughput is significantly low which makes it less 
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attractive for practical implementations. The PMF of the scheduled user location can then be given as: 



D. Location based Round Robin Scheduling Scheme 

Location based round robin is another non-opportunistic scheduling scheme which do not require any 
channel state information, however, it requires the location information of the users. Even though location 
based scheduling is not common in practice, the location of each mobile user can be determined at the 
BS using global positioning system (GPS) or estimate based on a power measurement of pilot signals 
from the surrounding beacons, e.g., using triangulation based techniques. In this regard, there are variety 
of techniques available in the literature which demonstrate how the location of users can be evaluated at 
the BS (see [[T3l . [[T4ll and the references therein). Moreover, the users located in different circular regions 
can also be classified based on the long term average SNRs, i.e., by computing SNR thresholds for 
different distances (rings) []4] which is a common technique in fractional frequency reuse (FFR) systems 
to distinguish between cell-edge and cell center users. 

In this context, we consider W time slots during which the distance of the users from their serving 
BS will remain approximately the same. For simplicity, the number of time slots W is set equal to K. 
At a given time slot T„,, we select any arbitrary user from a specific ring (analysis) and circular region 
(simulations) starting from the cell center. We continue to allocate the users by accessing the circular 
rings sequentially from cell center to cell edge. At this point, it is important to emphasize that all cells 
are considered to be time synchronous in allocating the users from particular rings, i.e., at a given time 
slot all cells are selecting the w*^ ring. Thus, the PMF of for a given time slot w denoted by 
P{r sei = can be given as: 



The scheme can produce relatively high capacity gains on average compared to the traditional round robin 
scheme. Moreover, the average fairness measure is very close to the traditional round robin. 



Uk 

U 



(16) 




E. Greedy Round Robin Scheduling Scheme 

Greedy round robin is an opportunistic scheduling scheme which captures the multiuser diversity while 
maintaining some fairness among users. In this scheme, we consider W = K time slots during which the 
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distance of the users from their serving BSs remain nearly the same, however, the small scale fading gain 
on the considered subcarrier may vary from one time slot to the other. We select the user with maximum 
SNR in each time slot Ty^, however once a user is selected from a ring, all users located in that ring will 
not be scheduled for transmission for the next K — 1 time slots. Note that all BSs are considered to be 
time synchronized in terms of scheduling. Clearly, the probability of allocating a ring k at Ti can simply 
be given by (fTTI) . However, the probability of selecting a ring k at T2 is a dependent event and can be 
derived using Bayes theorem as follows: 

P{r..^ = rl^) = j:P f 7^^^) ^(-sei = rf) (18) 

j=l \' sel — I j J 



P [ - r n (^c(-r..f ))"■ fc(7.^f ))"'" (19) 



where, 

p( 

Since the probability of allocating any ring k within time slot depends on all previous states, therefore, 
the principle of Markov chain transition probabilities is not directly applicable. For more clarity, the 
probability of selecting a ring k at T3 is given as follows: 

P(r,. -r^)-ttp ^ (^) nr.. ^ ^S) (20, 



3^k 



m 



where. 



In general, the probability of selecting any ring /c at a time slot T^, i.e., P(rsei = rl"") can be written as: 

y^y^... \^ p( = ^ ...pf ^sei = r/^ \ ^ ^ 



n^k s^n, j^s,n, 
sytk j^k 



(22) 

Computational Efficiency: The time complexity of the greedy round robin scheme is heavily based on 
the computational time of the N Integrate operation in MAT HEMATIC A. One N Integrate operation 
requires around 0.95 sec which is equivalent to the computational complexity of (i) greedy scheme and 
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(ii) the^r^f time slot of the greedy round robin scheme. Monte-Carlo simulation time required for 100,000 
trials in greedy scheme requires around 150.671928 sec which is more than 150 times the computational 
complexity of NIntegrate operation. This fact demonstrates the computational efficiency of greedy 
scheme in comparison to Monte-Carlo simulations. However, in the second time slot, greedy round robin 
scheme requires K — 1 integrations whereas, for the third time slot K — 2 integrations are required. 
Therefore, in general the computational complexity of greedy round robin scheme at any time slot w can 
be given as: 

Computational time[sec] = 0.95 + 0.95(i^ - z + 1); w <K (23) 

i=2 

where K denotes the number of rings. Therefore considering w = 15 and K = w for greedy round robin, 
the analytical time complexity is around 113 sec which is still lower than the Monte-Carlo simulation time 
required for the greedy scheme. Therefore, even though the greedy round robin scheme is computationally 
complex for large W, the evaluation time remains comparable to the Monte-Carlo simulations. 

F. Evaluating the Joint PMF of r^ei and 9 

Note that, P{rsci = rk) derived for all of the above scheduling schemes is the marginal PMF of 
-Pl'^sei = rk,0 = 6i) where 6 denotes the angle of the allocated user with respect to the serving BS 
and it is uniformly distributed from to 2n (see Fig. 1). Although the PDF of 6 is continuous we can 
discretize it for analytical consistency and complexity reduction. Consider discretizing the range of RV 
6^ in X uniform angular intervals of desired accuracy. Thus P{0 = Oj) = j where di denotes any discrete 
value that the RV 6 can take. Since rgd and 6 are independent, their joint PMF can be written as follows: 

P(rsei = rk,e = e;) = P(r,ei = rk)P{e = e,) = ^^""^y '''^ (24) 

IV. Part I: Distribution of Intercell Interference from One Cell 

The derivation for the distribution of the ICI from an interfering cell /, i.e., fx,{x), depends on the 
distribution of the distance between the allocated user in the interfering cell / and the BS of interest, i.e., 
ffseiji'^)- mentioned earlier, each interfering cell is assumed to have identical conditions in a given 
time slot. Therefore, ff^^^if) applies to all interfering cells and, thus, we will drop the subscript / in the 
sequel to simplify notation. Using the cosine law (see Fig. [U, we can write: 



^sei = rL + D'^- 2r,,iD cos9 



(25) 
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fsei is the distance of the allocated user in the interfering cell / from the BS of interest, rgci is the distance 
of the allocated user from its serving BS, i.e., (BS /), 9 E [0, 2n] and D = 2R since we consider universal 
frequency reuse with one tier of interfering cells. The approach can be extended to any number of tiers 
in a straightforward manner. In order to determine the PMF of fsd where fsd E [D — R, D + R], first of 
all we compute fi^k for given 6i and rk using (|65l) as follows: 

flk = rl + D^ -2rkD cos9i Wn,\/9, (26) 

where, fj ^ denotes the interfering distance from a specified polar coordinate (r^, 9i) in the interfering cell 
to the BS of interest located at a distance D (see Fig. [T]). In addition, it is worth to mention that fj ^ are 
the points at which P^fsd = fi^k) can be defined using (|63l) as follows: 

T->f~ ~ \ PyscX = ^k) /T7\ 

P[rsc\ = Ti^k) = ^ (27) 

The two dimensional data set of fsd, at which P{f sci = ?^i,fc) is defined, can then be grouped into M 
segments of any arbitrary width A. This can be done by dividing the distance between D — R and D + R 
into M equal segments of width A and mapping ft^k accordingly. Clearly, by adding all the probabilities 
for which fsd lies in the m}^ segment we get the probability of fsd = fm'- 

Pihel = fm)= ^ P{fsc\ = fi^k) (28) 

where fm denotes any discrete value that the RV fsd can take. Recall X = Kf~^x^ therefore the PDF of 
X conditioned on fsd can be determined by RV transformation as follows: 



jx\f,,i - — -—s — (2y) 



sel 

Averaging over the PMF of fsd, the distribution of the ICI, fx{x), from any cell / can be given as: 



fx{x) = V ^^^f^4-^P(fsel = fm) (30) 
- — - KVm 

It is important to emphasize that the derivation of the distribution of ICI is based on the scheduling 
decisions of interfering cells at a given time slot. Therefore, the parameter fm of the ICI distribution 
varies from one time slot to the other for the location based round robin and greedy round robin schemes. 
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V. Part I: MGF of the Cumulative Uplink and Downlink Intercell Interference 

Computing the distribution of the cumulative ICI Y requires the convolution of the PDF of L RVs Xi, 
V/ = 1, 2, ■ ■ ■ L, which is a tedious task for many practical scenarios. To avoid the convolutions, we utilize 
an MGF based approach and derive the expression for the MGF of the cumulative ICI Y. 

A. Derivation for the Uplink Interference 

Since each cell is considered to have same scheduling scheme deployed, therefore, the MGF of the 
cumulative interference considering i.i.d. interferers can be calculated as follows: 

L 

My{t) = X[MxXt) = {Mx{t))'' = (E[e*1)'' (31) 
1=1 

Looking at the structure of (1301 ), we can write Aix{t) as: 

/•OO p/~ f \ poo 

Mx{t)= e'^fx{x)dx = Mx{t) = Yl t~-P / e*Vx(^rti?-^)rfx (32) 

The derived expression is generic and applies to any composite fading distribution. Next, we will present 
explicit MGF expressions for the uplink ICI considering three typically used practical fading models. 
Special Case 1: Rayleigh fading -C,x ~ Exp(A): In this case, the small scale fading coefficient on a 
given subcarrier is considered to be Rayleigh distributed whereas the effect of shadowing is not considered. 
The distribution of interference considering a single interfering cell can then be derived as follows: 

fM 

fx{x) = K-'\rie-^''-^-'^P{hoi = r^) (33) 

f,n=fl 

Note that (1331 ) is a Hyper-Exponential distribution with parameter K~^\f^. Thus, using the MGF of the 
Hyper-Exponential distribution, Mvit) can be derived as follows: 

^- R^t'''"' - ) '''' 

Special Case 2: Generalized-/C composite fading ~ A^g("^s; "^c, In wireless channels, shad- 
owing and fading across the channel between a user and BS can be jointly modeled by a composite fading 
distribution. A closed form composite fading model, namely Generalized- /C also referred to as Gamma- 
Gamma distribution, has been recently introduced in [fT5l which is general enough to model well-known 
shadowing and fading distributions such as log-normal, Nakagami-m, Rayleigh etc. Using (|30l ), fxix) in 
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this case can be derived as follows: 

mc+ms 



fx{x) = J2 b^XriK-^ - Pinel = fm) (35) 



rm=ri 



where, K,;(.) denotes the modified Bessel function of second kind with order v, b = 2^ "'^^° . Performing 
some algebraic manipulations and using [16,, Eq. 6.643/3], the expression for Aixit) can be derived as: 

^ f fH^\ i>2,/5 /_ft2^/3\( 2 ) 



rm=ri 



where, W denotes the Whittaker function. Finally, A^y(t) can be written as follows: 



msj-2Tic--l\\ L 



(~ / ttiq +mr — J- \ 



Integrating the CDF and MGF of Generalized- /C RV which involves Meijer-G and Whittaker functions, 
respectively, in MATHEMATICA and MAPLE can be time consuming. 

Special Case 3: Gamma Composite Fading -C,,X ~ Gamma(ms, mc): In ifTSl . the authors proposed an 
accurate approximation of the Generalized- /C RV by the more tractable gamma distribution using moment 
matching method. The approximation provides a simplifying model for the composite fading in wireless 
communication systems. Using (l30l) . fx{x) can be written in this case as: 



rM 



Performing some algebraic manipulations and letting y = — t), Mx{t) can be derived as follows: 

(39) 

Finally, A^y(t) can be written as follows: 

H E - -J Gfi?£^ ,40, 
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B. Derivation for the Downlink Interference: 

The downlink interference Xi considering a single interfering cell / can be written as: 

Xi = f-\i (41) 

where xi is the interfering statistics from the l^^ neighboring BS, and f/ is the distance of /*'^ interfering 
BS from the scheduled mobile receiver. At this point, it is important to highlight that for a given and 
9i of a scheduled mobile user in the cell of interest, the distance of all interfering BSs can be calculated 
using cosine law which is not the same as in the uplink, where all fi are independent from one another. 
Therefore, conditioned on the location of the mobile receiver, the distribution of the downlink cumulative 
ICI, i.e., F|rfc, 6i, is simply a weighted sum of the distribution of the interfering channel statistics x- More 
precisely and by using the symmetry of the grid model the PDF of the cumulative ICI (Y = J2f=i -^i) 
can be given as: 

rK 2iT 

fy(y) = E E = ^i,k)fY\r„9Mrk, Oi) (42) 

rk=ri e,=0 

where 

L 

Y\n,ei = Y,ri''xi (43) 
1=1 

and the distance of all interferers can be determined using cosine law as follows: 



fi = ^D^ + rl- 2nDcos{{l - l)7r/3 + tt/Q - Oi) V/ = 1, ...L (44) 

Given the distance of the interferers, the conditional MGF of the cumulative interference can be calculated 
as follows: 

L 

MY\r,,e.{t) = \[M^{fih) (45) 
1=1 

Finally the MGF of the cumulative ICI can be calculated as follows: 

■^yit) =Y.Y1 = ^^,k)l[M^{f;^t) (46) 

rk=ri ei=0 1=1 



The expression for the MGF of the cumulative ICI in (1461) is general for any kind of composite channel 
fading models. The explicit expressions for three above discussed practical fading models can also be 
obtained in a straightforward manner. 
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VI. Part I: Evaluation of Important Network Performance Metrics 

In this section, we demonstrate the significance of the derived MGF expressions of the cumulative 
ICI in quantifying important network performance metrics such as the outage probability Pout, ergodic 
capacity C and average fairness among users numerically. 

Evaluation of Outage Probability: The outage probability is typically defined as the probability of 
the instantaneous interference-to-signal-ratio to exceed a certain threshold. In order to evaluate Pout 
numerically, we use the MGF based technique introduced in [|171 for interference limited systems. Firstly, 
we define a new RV, Z = qJ2f=i-^i ~ -^o = — ^o, where q is the outage threshold and Xq is 
the corresponding signal power of the scheduled user in the central cell. An outage event occurs when 
p{Z > 0), i.e., when the interference exceeds the corresponding signal power. This decision problem 
is solved in ifTTl by combining the characteristic function of Z and residue theorem. The characteristic 
function of Z is defined as (j)z{joj) = E,[e^^'^]. Considering interference Y and signal power Xq to be 
independent, the expression for (f)z{jui) can be given as, (f)z{jui) = 0y(jgu;)0Xo("~J'^); where (f)Y{qjuj) 
can be given by (l34l) . (l37l) . and (|40|) for different fading models. In general, the characteristic function of 
Xq can be calculated as follows IIT8I : 

roc roo 

(j)xS^) = E(e^^-^") = / e^'^^«/xo(a:o)rfxo = ju / e'''^'^ Fx,ixo)dxo (47) 

Jo Jo 

where Fxq{xo) = Y[f=i^iki^o) opportunistic schemes and compact closed form expressions of 
(f)Xo{juj) are available in the literature [|T9l Eq. 19]. For non-opportunistic scheduling schemes (pxoi^) = 
Y7r'^=ri^C\rkij^)'Pi^sc\ = ^k), whcrc ^^c characteristic function of C in ring k. The outage 

probability can then be computed by using the classical lemma introduced in [fTTl as follows: 

P„„..UirM(M^)^ (48) 

where lm(0^(a;)) denotes the imaginary part of (pzi^^)- Using (l48l) . the outage probability can be evaluated 
using any standard mathematical software packages such as MATHEMATICA. 

Evaluation of Ergodic Network Capacity: Another important performance evaluation parameter is the 
network ergodic capacity C, i.e., 

C = E log, I 1 + , ^° ^ 1 (49) 
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Usually, the computation of (|49l) requires (L + l)-fold numerical integrations. To avoid this, we utilize the 
efficient lemma derived in [|20l with a slight modification to take thermal noise into account and compute 
C as follows: 



E 



In 11+ ^" 



MY{t)-MxMt)^.-.H 







t 



e-" 'dt (50) 



where, A^y(t) = E[e-*^^=i^'] and Mxoxi^) = E[e"*(^o+Ef=i^i)] = E[e~*(^"+'^)]. Note that this is the 
definition of MGF as defined in [i20il which is not the same as our definition. Thus, we can use A^y(t) 
from (l34t . (|37] ) and (|40] ) directly with a sign change of jw. Moreover, (|49] ) can also be solved efficiently 
by expressing it in terms of the weights and abscissas of a Laguerre orthogonal polynomial [|20l : 



E 



In I 1 + ^' 



A Myi^.)-MxMQ , r, .... 
} 7 h-Ks (51) 



where and are the sample points and the weight factors of the Laguerre polynomial, tabulated in 
[fT6l . and Re is the remainder. Note that the MGF of Xq can be calculated as explained in (l47l) . 
Evaluation of Average Fairness: In order to quantify the degree of fairness among different scheduling 
schemes, we use the notion developed in ETIl . The average faimess of a scheduling scheme with U users 
can be given as, = — Yl'i=iPi ^iog°^u ' "^^ere Pi is the proportion of resources allocated to a user i or 
the access probability of user i. A system is strictly fair if each user has equal probability to access the 
channel and in such case the average faimess becomes one. The other extreme occurs when the channel 
access is dominated by a single user; in such case, the average faimess reduces to zero. The average 
faimess can be easily computed using our derived results as follows: 

where Uk denotes the number of users in a ring k. 

VIL Part I: Numerical and Simulation Results 

In this section, we first define the system parameters and describe the Monte-Carlo simulation setup 
which is required to demonstrate the accuracy of the derived expressions. We then address some important 
insights and study the performance trends of different scheduling schemes. 
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A. Parameter Settings and Simulation Setup: 

The radius R of the cell is set to 500m and the cell is decomposed into non-unifomi circular regions 
of width Afc. The path loss variation within each circular region is set io k = 2dB. For each Monte-Carlo 
simulation trial, we generate U uniformly distributed users in a circular cell of radius R. Each user has 
instantaneous SNR given by ([T]) and short term average SNR (7). We allocate a user with maximum 
instantaneous SNR in the greedy scheme whereas in the proportional fair scheme we allocate a user 
based on the maximum normalized SNR. For the round robin scheme, we select any user arbitrarily. For 
location based round robin we select a user randomly from the w*^ ring in a time slot w whereas we 
select a user with maximum SNR without considering the users of the previously allocated w — 1 rings 
for the greedy round robin scheme. Next, we calculate the distances of the selected users, i.e., r^ex, from 
their serving BS and compute the distances to the BS of interest, i.e., fsci for all scheduling schemes. The 
process repeats for large number of Monte-Carlo simulation trials. The distance data is then analyzed by 
creating a histogram whose bins are given by [1, ■ ■ ■ r^-i, r^, rfc+i, ■ ■ ■ R\. 

B. Results and Discussions: 

Fig. [2] depicts the PMF of the location of the scheduled user in a given cell based on the proportional 
fair, greedy and round robin scheduling schemes. Since the proportional fair scheme exhibits some fairness 
among users in a cell, the PMF of the allocated user locations is expected to be more flat compared to the 
greedy scheme. Since the cell edge has more users due to the large area and each user has equal probability 
to be allocated on a given subcarrier, therefore the round robin scheme exhibits high probability at the 
cell-edge. It is important to note that the numerical results for the derived PMF in Fig. |2] nearly coincide 
with the exhaustive Monte-Carlo simulation results with a small number of rings K = 10 and k, = 2dB. 
Moreover, it can also be noticed that the width of the circular regions tend to increase from cell center to 
cell edge which is due to the exponentially decaying path loss as mentioned in Section II. The number of 
required rings is expected to decrease by reducing (3 and increasing the amount of power decay within 
each circular region and vice versa. 

Another important point to explain with reference to Fig. |2]is that with the increase in the number of 
competing users on a given subcarrier, the PMF of opportunistic scheduling schemes tends to get skewed 
which is due to the fact that the higher the number of users in the cell center, the higher is the probability 
of allocating a subcarrier in the cell center. In order to get an integer number of users within a ring. 
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we perform rounding in the analysis, i.e., we consider zero active users in the rings where < 0.5. In 
Monte-Carlo simulations, we consider the probability of allocating a user in these rings to be zero which 
can also be verified from Fig. |2] 

In Fig. [3l the PMF of the distance between the allocated user in interfering cell / and the BS of interest, 
i.e., P(rsci = fm), is presented. Numerical results are found to be in close agreement with the Monte-Carlo 
simulation results. For the opportunistic scheduling schemes, it is likely that a user close to its serving BS 
can get a subcarrier, thus, the PMF of the distance of allocated interfering users is expected to have high 
density in the middle. However, the slight descend in the central region in Fig. |3]is due to ignoring users 
that lie within the rings where the average number of users is less than half. Moreover, we can observe 
that the round robin scheduler is highly vulnerable to interference compared to the other schemes as high 
interference is expected to come from the cell edge users in the interfering cells. On the other hand, the 
greedy scheduler is expected to have allocations near the cell center and, thus, leads to less interference 
from neighboring cells. The proportional fair scheme lies in between the two extremes. 

Fig. m illustrates the CDF of the ICI considering different number of interfering cells and path loss 
exponents (3 for the greedy scheduling scheme. With the increase in the number of interferers, the 
interference level increases. Moreover, as (3 increases, the signal degrades rapidly and thus interference 
level is reduced considerably. At this point, it is important to mention that in this paper we derive and 
utilize the MGF of the cumulative ICI rather than the CDF of the cumulative ICI in order to evaluate 
important network performance metrics. Therefore, the analytical part of the provided figure of the CDF of 
the cumulative ICI is plotted using a technique mentioned in ||22]| to convert MGF into CDF numerically. 

Fig. [5] investigates the effect of increasing the number of competing users on a given sub-carrier 
considering all scheduling schemes. It can be observed that the increase in the number of users enhances 
the performance of the opportunistic scheduling schemes due to additional multiuser diversity gains. 
The greedy scheme achieves the best performance whereas the round robin scheme achieves the worst 
performance. As expected, the proportional fair scheme lies in between the two extremes. The average 
capacity of location based round robin over W = K time slots has been shown to be better than 
the conventional round robin scheme. The average capacity results of the greedy round robin scheme 
is presented for W = 3 and W = 6. Clearly, for W = 1, the scheme is equivalent to the greedy scheme; 
however, with the increase of time slots, performance degradation takes place due to the reduction of 
multiuser diversity caused by ignoring the users from previously allocated rings. 
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Fig. [6] quantifies the average resource fairness of all presented scheduling schemes. As expected, round 
robin is a strictly fair scheme. The proportional fair scheme possesses the ability to enhance the network 
throughput compared to round robin scheduling while providing a high degree of fairness. The greedy 
scheme is observed to be the most unfair scheme. Considering K time slots, the average fairness of the 
location based round robin scheme is investigated and found to be very close to the round robin scheme, 
however, with degradation in performance as can be observed in Fig. |5l For the greedy round robin 
scheme, we plotted the fairness metric considering W = ?> and = 6; it is shown that as the number of 
time slots increases, the fairness improves with a trade-off price in terms of ergodic capacity. 

In Fig. |71 we evaluate the network outage probability as a function of the outage threshold; q = 
(Z + Xq)/Y for (i) U = 50 users; (ii) U = 100 users. The numerical and simulation results are nearly 
identical for most cases. The higher the outage threshold for a given signal and interference power, the 
greater outage is expected. Moreover, for larger number of users the outage probability is observed to 
reduce for all opportunistic scheduling schemes except the round robin scheme. Since increasing the 
number of users on a given subcarrier in non-opportunistic schemes does not directly affect the access 
probability of a ring k, therefore its impact on the ICI ia almost negligible. This fact can also be verified 
from Fig. |5] Finally, in Fig. [8l we evaluate the network ergodic capacity as a function of the fading 
severity parameter and average power of Gamma fading interference channels for different scheduling 
schemes. Firstly, it can be observed that increasing the average power of the interference channel which 
is given by = mcnis for a given fading severity parameter mg, the capacity degrades significantly for 
all schemes. Moreover, it is also shown that increasing the fading severity while keeping the average 
power = 3 fixed has minimal impact on the system capacity. Therefore, the lower average power of 
interference channel Vt, the better is the overall system performance. 

VIIL Part I: Conclusion 

We proposed a novel approach to model the uplink ICI considering various scheduling schemes and 
composite fading channel models. The proposed approach is not dependent on a particular shadowing and 
fading statistics, hence, extensions to different models is possible. The provided numerical results and 
help in gaining insights into the behavior of ICI considering different scheduling schemes and composite 
fading models. Moreover, they provide quantitative assessment of the relative performance of various 
scheduling schemes which is important for network design and assessment. 
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Part II of the Technical Report 

Part II: Introduction 

Energy efficient wireless communications has been gaining considerable attention these days mainly 
due to two major reasons i) dramatically varying global climate [|23l . and ii) slowly progressing battery 
technology Il24ll . In this context, power adaptation has been evolved as an efficient approach to reduce 
per capita power consumption, control inter-cell interference (ICI) and increase fairness among users in 
future generation wireless networks such as Orthogonal Frequency Division Multiple Access (OFDMA). In 
OFDMA networks, the system bandwidth is decomposed into orthogonal subcarriers. These subcarriers are 
adaptively allocated among users within a cell based on a predefined scheduling scheme and user transmit 
power levels. The allocated users on the same subcarrier in neighboring cells can cause significant uplink 
ICI depending on their transmit power level and channel conditions with respect to the base station (BS) 
of interest. 

Most of the recent literature considered the modeling of ICI in the downlink where the location of 
interferers is usually deterministic [|2], [|3]. However, compared to the downlink, the modeling of ICI in 
the uplink is more challenging due to the arbitrary locations of the interferers and the powers associated 
with them. Some interesting analytical models for uplink ICI have been presented in [?], [|6l; however, 
none of them considered the impact of channel based scheduling and power adaptation on the uplink ICI. 
Recently, in ll25l . we presented a semi-analytical framework to derive the distribution of uplink ICI on a 
given subcarrier assuming greedy scheduling without power adaptation. In this paper, we generalize the 
developed semi-analytical framework to incorporate the impact of power adaptation on the uplink ICI. 
This power adaptation promises considerable power savings while allowing high degree of fairness among 
users. Several power adaptation mechanisms are discussed in [|26l such as fast and slow power control, 
open-loop and closed-loop power control etc. In this paper we focus on slow power control considering 
that each user is capable of adapting its transmit power autonomously either by measuring its location 
through a global positioning system (GPS) or estimating its distance based on the power measurement of 
pilot signals from the surrounding BSs [[T4l . 

Part II: System Model 

We consider a given cell surrounded by L interfering cells. For analytical convenience, the cells are 
assumed to be circular with radius R. Each cell contains U uniformly distributed users where each user is 
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assumed to have perfect knowledge of its distance to the serving BS. The rate adaptation and allocation 
of users on a given subcarrier, therefore, depend on the channel qualities as well as the transmit powers 
of the users. The instantaneous SNR of any user can then be written as: 

min(Pmax, Por^)r-^C ,c 
7 = 2 ^^^^ 

where Pmax[W] is the maximum transmit power of a user, Po[W] is the desired power level at the receiver, 
r[m] is the user distance from its serving BS, /3 is the path loss exponent, cr^ denotes the thermal noise 
at the receiver which is considered to be unity without loss of generality and represents the combined 
shadowing and fading random variable (RV). More explicitly, (l53l) can be re- written as: 




7 = < (54) 



The distance at which users need their maximum power to compensate path loss completely is referred 
to as threshold distance (r^) and can be computed as follows: 



P. 



■^■nifj '''' 

Users located within can compensate path loss completely while saving some proportion of their power, 
whereas the users located beyond rt transmit with their maximum power to achieve a certain rate that is 
less than their desired target. 

Each cell is decomposed into K concentric circular rings. The circular regions between two adjacent 
rings are characterized by uniform path loss variation (in dB) and, thus, possess non-uniform width A^. 
Since path loss varies exponentially with distance, increases from cell center to cell edge. Thus, the 
number of circular regions in each cell depends on the path loss exponent. The average number of users 
in a given ring k can be computed as follows: 

^2 = 1,2, ■ ■ ■ ,K, (56) 

where denotes the radius of ring k. It is important to note that Uk can be a fraction of a number; 
therefore, we round off the fractional part of users in each ring. The motivation behind dividing each 
cell into a number of circular regions is that in each region the channel conditions of the users become 
relatively similar especially for large values of K. 
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The proposed approach to model ICI is detailed in the following steps: 

• Derive the distribution fr^^^ir) of allocating a given subcarrier to a user at a distance r^ei from its 
BS. 

• Derive the distribution of the distance between the allocated interfering users and the BS of the cell 
of interest, i.e., determine /r^^i(^) using fr^^^ir) where f^ei is the distance between interfering users 
and the BS of interest. 

• Derive the distribution of the ICI fxXx) from the allocated user in neighboring cell / to the BS of 
interest. Since the allocated interfering user can transmit with different power levels depending on 
the distance from its own serving BS, the incurred interference can be modeled as 

{PQr^f~^X fE[D—rt D+rt\,rE[0 rA 
(57) 
Pmaxf~^X otherwise 

where D = 2R and x denotes the combined shadowing and fading component of the interference 
statistics. 

• Derive the MGF of the cumulative interference Y = "^f^i Xi caused by the allocated interfering 
users in all neighboring cells. 

Part II: PMF of Allocated User Locations 

In this section, we derive the discrete distribution of the distance of the allocated users in a given cell, 
i.e., the probability mass function (PMF) of Vsci- The derivation is divided into two steps explained as 
follows: 

Step 1 (Selecting the user with the highest SNR in ring k): Since each circular region has uniform path 
loss variation, the users within a ring k are assumed to be subject to approximately the same path loss. 
Thus, selecting a user in a ring k is equivalent to selecting the user with maximum channel gain among 
all the users in ring k, i.e., 

Cfc = max{Ci, C2, ■ ■ ■ , 0, ■ ■ ■ , C«J (58) 

For simplicity, we consider independent and identically distributed (i.i.d.) channel gains of all users. 
Therefore, for any ring k, the CDF and PDF of the maximum channel gain Ck can be written as follows, 
respectively: 

F,ACk) = \{F,^{Ck) = {F,{Ck)r (59) 

i=l 
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Considering the model in (|54] |. we split the analysis into two regions, namely the region within the 
threshold distance and the region beyond the threshold distance. After performing the RV transformation, 
we can write the CDF of the selected user SNR in each ring k (7^) as follows: 

^7.(7.) =r ' (61) 



Step 2 {Selecting the ring k with maximum SNR from the K rings): In this step, we compute the probability 
of selecting k^^ ring among all other rings. It is important to note that this is equivalent to selecting the 
ring k which possesses the user with the highest SNR among all rings. Thus, conditioning on 7^., the 
PMF of Tsci can be written as follows: 

K 

P(rsei = rfc|7fc) = J]F^,(7fc) (62) 

i=l 

By averaging over 7^, the final expression for the PMF of Tsci can be written as follows: 

/"OO 

^('^sei = Tk) = / P(rsci = rfe|7fc)/^^(7fc)rf7fc (63) 
Jo 

The result in (|63] ) can be evaluated accurately using standard mathematical software packages such as 
MAPLE and MATHEMATICA and is valid for any composite fading statistics. 

Note that P(rsci = r^) in (|63l) is the marginal PMF of P(rsei = rk,9 = On) where 6 is the angular 
position of the allocated user and is distributed uniformly from to 2tx. Although the RV 6 possesses 
continuous distribution, we discretize it in order to reduce complexity. Thus, discretizing 6 into N uniform 
angular intervals, P{6 = On) is where 6'„ denotes any discrete value that the RV 6 can take. Since 
Tsex and 6 are independent, their joint PMF can be written as: 

P(rsci = rk.e = On) = (64) 

Part II: Distribution and MGF of the ICI 

In this section, firstly we find the distribution of the distance of the users allocated in the interfering 
cell / to the BS of interest, i.e., ff^^^ (f). Based on the derived expression we derive the distribution of 
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ICI from / interfering cell, i.e., fxi{x) and the MGF of the cumulative ICI Y. 

A. Distribution of the allocated interfering user locations 

Since, each cell is assumed to have identical conditions, ff^^i{f) remains the same for all interfering 
cells and we will not use subscript / any further to simplify notations. Using the cosine law, we can write: 

^sei = ^'ei + D^- 2r,,iD cose (65) 

where fsci is the distance of the selected interfering user in cell / from the BS of interest, r^ci is the 
distance of the selected interfering user from its own BS, i.e., (BS /), 9 E {0, 2n}. 

In order to determine the PMF of fsd where fgci G {D — R, D + R}, first of all we define f„ ^ for given 
9n and using (1651 ) as follows: 



fn,k{r, 0) = ^Jrl + - 2rkD cosOn Vr^, Ve„ (66) 
Clearly fn,k{r, 0) are the points at which P{f sci = ?^n,fc) can be defined using (|63l) as 

P{rsci = rn,k> = (67) 

B. Part II: Distribution of the ICI from one interfering cell 

Since the interfering users can transmit with different power levels depending on their distance from 
their serving BS, the interference X can be categorized into two regions mentioned as follows: 

X={ ' (68) 

^max^~fcX otherwise 

where x denotes the interference channel statistics. The PDF of X conditioned on fn,k{r,9) can be 
determined by RV transformation as follows: 



r. 



n,k e[D-rt D + rt],rk E [0 n] 



fxKJx) = { _^ ^ " " (69) 



^fxip^'^nk) otherwise 
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Simply averaging over f„^fc and letting A = ■^r^ki''^,0)r^'^ and B = ■p^r^ki''^,^) we can write the 
distribution of interference, i.e., fxi^) as shown below: 

, , , , Er-, Ee„ ^ ^(^sei = fn,k)fxiAx) G [Z^ - ^ + r^], r, G [0 

= { (70) 

Erfc Ee„ BP{f^^x = fn,k)fx{Bx) otherwise 
Finally, fx{x) can be written explicitly as follows: 

fx{x)= 5Z^P(rsei = r„,fe)/^(Ax)+ $^5P(rsci = r„,fe)/;,(i?x) (71) 

MGF of the cumulative ICI 

Computing the distribution of the cumulative ICI Y requires the convolution of the PDF of L RVs Xu 
V/ = 1, 2, ■ ■ - L, which is a tedious task for many practical scenarios. To avoid the convolution operations, 
we utilize an MGF based approach. Since the scheduling scheme is considered to be identical in all cells, 
the interferers are i.i.d. and therefore the MGF of the cumulative interference Y can be written as follows: 

L 



Mrit) = l[MxXt) = {Mxit))'' = (E[e*^])'' (72) 



1=1 

Looking at the structure of dTOl ). we can derive MGF of any composite fading model as Aixit) = 
/o°° e^'-^fx{x)dx The expression applies to any kind of composite fading models. Due to space limitations 
we will study only the MGF of the Gamma composite fading case, i.e., we consider a scenario in which 
shadowing and fading statistics are modeled by a Gamma and Nakagami distribution (also referred as 
Generalized- K IfTSlI ). respectively. Recently, in IfTSl an accurate approximation of the Generalized- K 
RV using moment matching method has been proposed to increase its analytical tractability, i.e., the 
Generalized-K distribution can be approximated by a simple Gamma distribution [[TSl . Therefore in this 
case Aix{t)) can be derived as follows: 



(73) 



fk&lo n] 9 ^ " ^ r^([0 rt] 0„ ' 

Finally A^y(t) can be given simply using (|72l) . 

Part II: Evaluation of Network Performance Metrics 

In this section, we will utilize the derived expressions to evaluate the network ergodic capacity, average 
faimess, and average power preservation per user. The evaluation of outage probability is skipped due to 
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space limitations, however, the readers can refer to [|25l for details. 

C. Evaluation of Network Ergodic Capacity 

Using the lemma derived in [|20]| . the following expression is valid for interference limited systems: 



E 



In 1 + 



My[t)-M.Mt),^ (74) 
t 



where, Myit) = E[e'*^^=i is the MGF of the cumulative interference and Mxo,Y{t) = E[e"*(^"+^^=i 
]g|^g-t(Xo+y)j |g jj^g joint MGF of the corresponding signal power of the scheduled user Xq and cumulative 
interference Y. Since Xq and Y are independent, A^Xo,y(^) = A^j)s:o(^)-A^y (^)- The expression for Aixoi^) 
can be given as: 

POO POO 

Mxo{t)= / e'""' fxoixo)dxo = / te*"°Fxo(xo)rfxo (75) 
Jo Jo 

where Fx„(a;o) = Y[f=i ^-/ki^o), fxai^o) = '^^Xoi^o)- Closed form expressions are also available for 
Mxoit) in the literature llT9l . 

D. Evaluation of Average Fairness and Power Preservation 

To measure the degree of fairness among users, we follow the notion developed in [|27l . The average 
fairness in a network with U users is defined as = — Yl'i=i Pi^^j where Pi is the proportion of resources 
allocated to a user i or the access probability of user i. A system is strictly fair if each user has equal 
probability to access the channel and in such case the average fairness becomes one. The average fairness 
can be easily computed using our derived results as: 

= - V P{rsei = rk) — (76) 

t^i logf/ 

where P(rsci = r^) is given by (l63l) . Moreover, the average power savings per subcarrier can also be 
calculated as follows: 

n 

P=Y1 ^(^-1 = ''k) (^max - Porf ) (77) 

where rj denotes the threshold distance. 

Part II: Results and Analysis 

In this section, we aim to validate the accuracy of the derived expressions through Monte-Carlo 
simulations. The results are presented for Gamma composite fading, i.e., C ~ Gamma(l, 1), x ~ 
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Gamma(3/2, 2/3). Firstly we will provide a brief overview of the Monte-Carlo simulation setup. 

1) Generate U uniformly distributed users per cell. Each user has an instantaneous SNR given by (|53l) . 
Select a user with maximum instantaneous SNR. Store the distance of the selected user, i.e., Vsci 
from the serving BS. 

2) Compute the distance of the selected user from the BS of interest, i.e., f^.^ using cosine law and 
finally generate the interference using (1571 ). 

3) Repeat all steps for a large number of iterations. Generate histogram for the discrete RV rsci with 
non-uniform bin widths. 

In Fig. |9] the impact of the maximum transmit power limit is shown on the PMF of allocated user 
locations. The obtained PMF results fit well with exhaustive Monte-Carlo simulations. Since the users 
located within the threshold region can compensate their distance based path loss, each user has on 
average equal probability of allocation within r^. The increasing trend of PMF within rt is therefore 
simply due to an increase in the number of users in each ring from cell center to the cell edge. It is 
important to note that the users located beyond transmit with their maximum power as they cannot 
compensate path loss. These users are therefore scheduled based on their relative channel gains which 
prioritizes close users over the far users and hence causes rapid decay of allocation probability beyond 
Tt- In greedy scheduling [|25l . the cell center users have higher priority to be allocated over the cell edge 
users. On the other side, round robin scheduling provides equal probability of allocation to each user, 
hence high probability of allocation near the cell edge due to the large area and large number of users 
at the cell-edge. By observing the result, it can be concluded easily that greedy scheduling with power 
control (PC) follows the trend of round robin within rt whereas the trend of greedy scheduling beyond r^. 
The performance of greedy scheduling with PC is therefore expected to lie in between the two extremes. 

Two different transmit power limits are also studied in Fig. [9] which yields two threshold distances, 
i.e., rt = 400m and rt = 260m, respectively. It can be observed that the greater the maximum transmit 
power, the greater is the threshold distance and more users located farther from the BS become capable 
to compensate path loss which increases fairness and in turn the incurred ICI. The slight mismatch in the 
simulations and analysis demonstrates the impact of discretization which is dominant for channel based 
scheduling beyond r^. However, this error can be reduced by increasing the number of rings. 

The CDF of the ICI for different transmit power budgets and different path loss exponents for greedy 
scheduling with and without PC is plotted in Fig. [TOl High values of path loss exponents causes rapid 
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signal degradation, hence, reduces ICI. Moreover, it can be observed clearly that with low user transmit 
powers, there is a significant reduction in ICI compared to the high transmission powers. It is further 
interesting to note that the performance of greedy scheduling with PC always remain better than the greedy 
scheme in terms of incurred ICI, average fairness (see Fig. WVpartia)) and average power consumption of 
the users. The top figure in Fig. [TT] quantifies the average fairness of the greedy with and without PC and 
round robin schedulers. With the increase of transmit powers, far users can also adapt their power which 
increases the average fairness among users. For high user transmit powers, the greedy scheduling with 
PC achieves the fairness of round robin scheme as is also evident from Fig. |9] However, it is important 
to note that the capacity (see Fig. [TT]3art(b)) and power preservation remains better than the round robin 
scheme in which power savings are zero. 

The bottom figure in Fig. [TT] demonstrate the network capacity of interference limited systems (i.e., 
thermal noise is neglected). Without PC, the performance of greedy and round robin scheduling remains 
independent of the transmit power as the factor of Pmax cancels out in the capacity calculation. However, 
since the greedy scheduling with PC have less ICI then the greedy scheduler, the network capacity is 
expected to increase which is not the case as the corresponding user transmit powers are also lowered along 
with the interfering powers. The main reason of the capacity degradation with the increase in transmit 
power budget is that the greater transmission power more users can compensate path loss which reduces 
the number of users transmitting with their maximum powers. This phenomena on one hand increase 
average power savings whereas on the other hand degrades overall system capacity. 
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Fig. 1. Geometrical illustration of dividing the cellular network into multiple rings of non-uniform width Ak- 
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Fig. 2. PMF of the distance of the allocated users in a given cell (i.e., PMF of rsci) for proportional fair, greedy and round robin scheduling 
schemes with path loss exponent P = 2.6, U = 50, Number of Monte-Carlo simulations =100,000, C=60 dB, Pniax=lW, a^=-174 dBm/Hz. 



34 



"0.06 
0.04 
0.02 



Proportional Fair Scheduling 



Histogram by Simulations 



- • - Analytical PMF 




0.2 



0.1 



§ 

■f= 600 700 800 900 1000 1100 1200 1300 1400 1500 

c Greedy Scheduling 

1^0.2 

CO 
CO 
CC 

^0.1 



■ -■- 



■ ■ — ■ — ft — * 



Id 

5 600 700 800 900 1000 1100 1200 1300 1400 1500 

o 



Round Robin Scheduling 





600 700 800 900 1000 1100 1200 1300 1400 1500 

Distance of the allocated interfering users from the BS of interest [m] 



Fig. 3. PMF of the distance at which the users in the interfering cells are allocated (i.e., PMF of fsei) for proportional fair, greedy and 
round robin scheduling schemes with path loss exponent /3 = 2.6, U — 50, T = 720, x ~ Gamma(3/2, 2/3), Number of Monte-Carlo 
simulations =100,000, C=60 dB, P^^^=m, (j'^=-\lA dBm/Hz, A=50 m. 
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Fig. 4. Impact of different number of interferers L and various path loss exponents (/3) on the CDF of cumulative ICI considering greedy 
scheduling scheme, U = 50, T = 720, x ~ Gamma(3/2, 2/3), Number of Monte-Carlo simulations =100,000, C=60 dB, P^^^=1W, 
a^=-n4 dBm/Hz. 
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Fig. 5. Network ergodic capacity for different number of users considering different scheduling schemes with /3 — 2.6, x 
Gamma(3/2, 2/3),C=60 dB, P,^^^=m, a'^=-\lA dBm/Hz. 
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Fig. 6. Average system fairness for different number of users considering different scheduling schemes with /3 = 2.6,C=60 dB, Pmax=lW, 
(7^=-174 dBm/Hz. 
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Fig. 7. Impact of different scheduling schemes on the network outage probability for different number of users and various outage thresholds, 
13=2.6, X ~ Gamma(3/2, 2/3), C=60 dB, P,^^^=m, cr^=-174 dBm/Hz. 
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Fig. 8. Ergodic network capacity of different scheduling schemes as a function of different parameters of the interference statistics 
X Gamma(ms,mc) with U — 50 users, 13=2.6, a'^=-m dBm/Hz. 
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Fig. 9. PMF of the distance of allocated users from their serving BS in a given cell (i.e., PMF of r^^\) for greedy scheduler with power 
adaptation, _R=500 m, /3=2.2, f/=50, Po=-23 dBm. 
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Fig. 10. (a) CDF of the ICI for different transmit power levels and path loss exponents considering greedy scheduling with and without 
power control (PC), Po = -23 dBm, R = 500 m. 
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Fig. 11. (a) Average fairness among users considering greedy scheduler with and without power control (PC) and round robin scheduler 
without PC (b) Ergodic capacity considering greedy scheduler with and without PC and round robin scheduler without PC, U = 50, Po = 
-23 dBm, i? = 500 m, /3 = 2.2. 



